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We compute the graviton induced corrections to Maxwell's equations in the one-loop and weak 
field approximations. The corrected equations are analogous to the classical equations in anisotropic 
and inhomogeneous media. We analyze in particular the corrections to the dispersion relations. 
When the wavelength of the electromagnetic field is much smaller than a typical length scale of 
the graviton two-point function, the speed of light depends on the direction of propagation and 
on the polarisation of the radiation. In the opposite case, the speed of light may also depend on 
■ the energy of the electromagnetic radiation. We study in detail wave propagation in two special 

' backgrounds, flat Robertson- Walker and static, spherically symmetric spacetimes. In the case of 

£sQ , a flat Robertson- Walker gravitational background we find that the corrected electromagnetic field 

! i ■ equations correspond to an isotropic medium with a time-dependent effective refractive index. For a 

O ' static, spherically symmetric background the graviton fluctuations induce a vacuum structure which 

. causes birefringence in the propagation of light. 

I/-) . PACS number(s): 04.60.-m, 11.15.Kc. 

I. INTRODUCTION 

> ■ 

ov . 

It is well known that when the QED vacuum is modified by external conditions such as background fields, finite 
£NJ . temperature, or boundary conditions, the propagation of the photons can be affected in a non-trivial way. The vacuum 
behaves as a dispersive medium in which the propagation of light generally depends on the direction of propagation 
and on the polarisation of the radiation. Physically, the effect can be understood as follows: the photon exists for part 
of the time as a virtual e~e + pair, on which the external conditions do act and modify the propagation. In previous 
works the problem of photon propagation in modified QED vacua has been analyzed for external electromagnetic 
(EM) fields fj]||, boundary conditions fj|,p|, external gravitational fields ||[|, an< ^ finite temperature 0. Further 
references and details can be found in Ref. [pj . There is also an experiment under construction to detect birefringence 
Oh, of the QED vacuum in the presence of a strong magnetic field M . 

The phenomenon is of course quite general, and not restricted to e~e + pairs. The interaction of the photon with 
any other field will produce similar effects: the photon will not follow, in general, a geodesic of spacetime []. In this 
paper we will analyze the effect of the coupling between a classical electromagnetic field and a quantum gravitational 
field on the propagation of EM radiation waves. We will show that, indeed, the graviton loop leads to effects that are 
similar to those already studied and calculated for QED vacua. 

Even though General Relativity is a non-renormalizable theory, the one-loop corrections are meaningful when 
the quantized gravitational field theory is looked upon as an effective field theory pL It is possible to compute, for 
instance, the leading (long distance) quantum corrections to the Newtonian potential |12| , |l"3| . Our calculation provides 
another example of a quantum gravity effect that can be estimated using General Relativity as a low energy effective 
field theory for quantum gravity. Moreover, it could also be of some interest from a phenomenological point of view. 
Indeed, Amelino-Camelia et al ]l4| pointed out that many quantum gravity scenarios predict a frequency-dependent 
velocity of light that could be observable for (cosmological) gamma-ray bursts. Gambini and Pullin jl5j studied the 
propagation of light in canonical quantum gravity and found that the modified Maxwell's equations imply a frequency 
and hclicity dependent velocity of propagation. We will see that, in principle, similar results can be found by taking 
into account the interaction between gravitons and EM radiation in the low energy theory. 
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1 It is even possible to have superluminal propagation. However, as extensively discussed in the literature ]^fl|lo|], this does 
not imply causality violations. 
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The paper is organized as follows. In Section II we obtain the quantum corrections to the classical Maxwell's 
equations induced by the graviton loop. We write the corrections in terms of the coincidence limit of the graviton 
two-point function. In Section III we present a qualitative analysis of the quantum corrections to the dispersion 
relations. We analyze the cases A >> L c and A << L c , where A is the wavelength of the classical electromagnetic 
radiation and L c is a typical scale of variation of the graviton two-point function. We show that in general the velocity 
of light depends on the direction of propagation and on the polarisation, i.e., we find gravitational birefringence. In 
Section IV we discuss quantitatively two examples: a flat Robertson- Walker (RW) background metric and a static 
metric with spherical symmetry. Section V contains our conclusions and final remarks. 



II. EFFECTIVE EQUATIONS OF MOTION 



Consider pure gravity described by the Einstein-Hilbert (classical) action f\ 



S G = 



A A x^/—gR 



where k 2 — 32irG and R is the Ricci scalar. The classical action for the EM field is given by 



(1) 



(2) 



where F„ 



V ' V A^ is the field strength tensor and the gauge potential. The classical energy-momentum 



tensor associated to the EM field is given by 



J EM 



(3) 



The classical action of the EM field depends on the (classical) gravitational background. It is then natural to ask 
ourselves what are the effects on Maxwell's equations due to a change in the gravitational background, and in this 
paper particularly, what are the implications of the one- loop graviton fluctuations. 

We define the classical action of the combined system (gravitational field plus classical EM radiation) 



s, 



clas 
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EM 



■s, 



G ■ 



(4) 



The effect of quantum metric fluctuations can be analyzed by means of the background field method, expanding the 
total action S^ias around a background metric as t^ t „ — > g^+Kh^, and integrating over the graviton field {h^ v ) degrees 
of freedom to get an effective action for the background fields g^ v and . In order to fix the gauge one must choose a 
gauge-breaking term x M [ff> h]> with its corresponding gauge-breaking action S gaugc [g, h] = —(1/2) J d A x^/—gx tl g^vX u , 
and its corresponding ghost action Sghost Pql - The complete effective action 5* ff is obtained by integrating the full 
action S = S c i aB + S gaugc + S g h os t over the graviton and ghost fields. We evaluate this effective action in the one-loop 
approximation, for which we expand S up to second order in gravitons. The second order term can be shown to have 
the form 17 



S (2) = J A A x^jh^{0^ aT +P^ aT )h 



(5) 



where O = 0^ uaT is a second order differential operator that depends on the background metric and is independent 
of the EM field |17j (we will not need its exact form in what follows), and the tensor P = p^ r7T arises from the 
expansion of S'em to second order in gravitons, and reads 



FopF"* [ -y v g° T - ^g VT 



-g !icr F OLT 



^F al 3g aT {~g^F av + g a ^F^) + 2F af} (F aT g fJ > i g"< T + F vfi g aa g^ T + F VT g a " g 0a ) 



(6) 



2 Our metric has signature ( h ++) and the Riemann and Ricci tensors, and the scalar curvature are defined as R^ ua( . 

daT^g — ■ ■ ., R a /3 = R^ a u,B' &n< ^ R = 9 al3 RaB, respectively. We use units such that h — c = 1. 



There is also a second order term in the ghost fields that for gauge-breaking terms linear in the metric fluctuations 
decouples from the gravitons and is only coupled to the background metric. This means that the one-loop effective 
action for the combined system reads [jl7| 

S e s[g^ ,Aft] = 5^ + -Trln(0 + P) - iTr In Gghost , (J) 

where Gghost is the second order differential operator that comes from integrating over the ghost fields. 

It is extremely complicated, in general, to calculate the one-loop effective action. In this paper we will use the weak 
field approximation, assuming that the EM field is a test field that does not affect the background metric g^. In this 
approximation the effective action takes the form 

S'cff = S clas + ^TrlnO - zTrlnG ghost + J &Xy/=gF ya ,(x)F aT {x)M'"' ir (x, x) , (8) 



where 



M* J,vcrT (x,x) = - — (h a 0(x)h Cn (x) + hc v {x)h a p(x)) 



+ lg aP (g K g™g Tri - s M - g ail g K g TV + sT'g^g"*) 
+ s a WV^ + g° g TV g m - W 3 - 9° W 7 ) 

+ <? CT V W C - <f W Q <7 K ] ■ (9) 

Note that the tensor M^ vaT has the following symmetry properties, which are similar to the symmetries of the 
Riemann tensor: M>" <7T = -M ufl(7T , M^ aT = -M* VTa , and M fIl " 7T = M OT ^ v . In view of these properties, the 
only non- vanishing components of this tensor are M 010 ^, M 0l: > k , and M y , where i,j,k,l are spatial indices. This 
tensor has 21 independent components and depends on the two-point function of gravitons, evaluated in an arbitrary 
quantum state |\&). The two-point function of gravitons can be written as follows 

G^ aT {x,x') = (^\h^{x)h aT {x') + h ar {x')h' ll '(x)\^) , (10) 

taken at coincidence (x = x'). In the following we will assume that the graviton state |\&) preserves the symmetries 
of the background metric, so that the tensor M*"" 71 will share those same symmetries. 

The (one-loop) gravitationally modified equations of motion for the EM field, SSgg/SA^ = 0, are given by 

V^G^ = , = F"" - AM^ aT F' TT . (11) 

These equations are analogous to the classical Maxwell's equations in the presence of a linear, anisotropic, and non- 
homogeneous media. To see this explicitly, we recall that in a local Lorentz frame F 0m = E m and F mn — e mn k B k , so 
that we can introduce the vectors G 0m = D m and G mn = e mn k H k , namely 

D* ee & - 8M% m E m - 4M\ n e m %B k , (12a) 



Om 

1 

The quantum corrected equations in a local Lorentz frame read 



H* ee B> - 2e jmn M mnik e ik Q B« - Ae jki M ki0m E m . (12b) 



V • D = , (13a) 

V • B = , (13b) 
VxH = 9,D, (13c) 
VxE = -d t B , (13d) 

where we have also included the equations that follow from the Bianchi identity V M F, CT + ^uF afJ _ + V aF^v = 0. 
The constitutive relations D[E,B] and H[E, B] are unusual when Moyfe ^ 0, since in this case both the electric and 
magnetic fields appear in the definition of D and H. However, when Moijk = the relations are exactly equivalent 
to those of a linear medium characterized by spacetime dependent electric permitivity and magnetic permeability 
tensors defined as 
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D l = e^Ei , e l 3 = S l } - 8M iQ 0j , (14a) 
B i = itjH* , ^ = 5 l 3 + 2e l mn M mnab e abj . (14b) 

We emphasize two relevant technical points. On the one hand, the graviton two-point function will depend, in 
general, on the choice of the gauge-breaking term for the graviton fluctuations. However, the background metric 
will also depend on the gauge-breaking term through the semi-classical Einstein equation. Both dependencies should 
cancel out in the dynamics of any test field (see Refs. jl3| and p8[). We will not consider this problem in what follows, 
i.e., we will assume that the background metric already contains the back reaction of gravitons, computed with a 
given gauge-breaking term. On the other hand, the two-point function will diverge at the coincidence limit. Adequate 
counterterms are needed to absorb the divergences. In the spirit of effective field theories, these counterterms will not 
contribute in the long distance/low energy limit, which will be dominated by the non-local, non-analytic part of the 
two-point function pi Jl3 1 . 



III. QUALITATIVE ANALYSIS 

Given the quantum corrected equations, we can distinguish two different physical regimes depending on the relative 
size of the wavelength A of an EM radiation field described by F pv and the typical scale of variation of M pvaT , L c . 

When A <C L c , we can take the tensor M pvaT outside the covariant derivative in Eq. (|ll]), as it does not change 
significantly over the scale A. The equation of motion can then be written in the form 

V M -F^ - 4M^ T V^F <TT = . (15) 

We first introduce the following variables ^ 

= f^e'* , (16) 

with the amplitude and ip the phase, such that fc M = V M V is the momentum of the EM wave. We assume that 
the amplitude f^ v is the slow varying variable and ip is the fast varying variable, so that from now on we discard any 
gradients and/or time variations of the amplitude / M „. 

We start from equation ( |l5| ) and make use of the new variables to write 

k p r - ^ p M pv ° T ] aT = . (17) 
The remaining Maxwell's equation implies the following 

kfifua + k u f afl + k a f pl , = . (18) 

We now multiply equation (|l7]) by A; M and make use of equation ( |l8| ) to obtain 

fc 2 /^ + Sk^M^Pk^frp = . (19) 

This equation is similar to the one discussed by Drummond and Hathrell in Ref. ||. (See also Ref. ||.) In those 
references the corrections are due to fermion loops and are proportional to the Riemann tensor R^ u P a (this is true in 
the case of empty spacetimes, so that R^" and R vanish ||), instead of M A " yp<T . 

In the absence of quantum corrections one obtains the usual dispersion relation k 2 = 0. The graviton loop induces 
modifications to this relation, i.e., light rays do not follow null spacetime geodesies. We will analyse in detail several 
examples in the following sections. Here we discuss qualitatively some general properties of the modified dispersion 
relations. 

It is easy to see that the tensor M^ pa is dimensionless and proportional to the square of the Planck length Lp ~ n. 
Therefore we expect M^ vp<T = (Lp/I/ c ) 2 C MI/pcr , with C^ 1 " 7 a dimensionless, slowly varying tensor. Inserting this in 
Eq. (|l9|), we see that the modified dispersion relation will have the general form 

k 2 + c^fc" = , (20) 

where c M „ is a slowly varying tensor of order 0(Lp/L 2 ) that depends on the direction of propagation and the 
polarisation of the EM radiation. Therefore, we expect the modifications in the speed of light to be proportional to 
Lp/ L 2 . and gravitational birefringence of the same order of magnitude. 

We consider now the opposite case A 3> L c . For simplicity, and in order to be able to compare with previous works, 
we assume that the background metric is flat and that the quantum state for gravitons is such that the two-point 
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function has a random variation on micro-scales (much smaller than any other scale of the system, but still larger 
than Lp). In other words, we assume that the spacetime looks classical at scales larger than L c and has a complicated 
random structure at scales smaller than L c . This kind of states has been considered before in the context of loop 
quantum gravity [|l5Ul9| . 

In this situation, the quantum corrected equation ( |Tl] ) reads 

- 4M^ T (d„F° T ) - 4(d^M^ aT )F° T = . (21) 

It is, of course, not possible to neglect the derivatives of the tensor M^ vpa . The equation describes the propagation 
of a classical electromagnetic wave in a random media. 

As the wavelength is much larger than L c , in order to obtain a modified dispersion relation we average the field 
equation over a spacetime domain of size L 4 , with L c << L << A. To compute the average of the products (dM)F 
and M{dF) we expand F around a point xq in the domain. Schematically F(x) « F(xq) + dF(xo)(x — xq). Denoting 
by (...) the average over the domain and using that for a random structure (M) — and (DM) — 0, we obtain 
(FdM) w dF(x )(dM(x — x )) and (MdF) w d 2 F(xo)(M(x — xq)). Therefore, in this approximation, the average 
of Eq. ( p?i"| ) will contain higher derivatives of the EM fields (as long as (M(x — xq)} ^ Q). As a consequence, on 
dimensional grounds we expect a modified dispersion relation of the form 

k -\- c^isk^k -(- c^ypk^k = , (2^) 

where c M „ = 0(Lp/L 2 ) and c^ p = 0(L 2 P /L). 

The quadratic correction modifies the speed of light as in the previous case. The cubic term is qualitatively different, 
since it produces a variation of the speed of light that increases linearly with the energy of the EM radiation. As 
already mentioned, this kind of correction may be relevant from a phcnomcnological point of view, because it induces 
a non-trivial structure in the arrival time of light rays coming from gamma-ray bursts Jl4|,|l5| . 

IV. EXAMPLES 

In the following sections we will concentrate on two particular classes of gravitational backgrounds, namely flat RW 
metrics and static spherically symmetric backgrounds. 

A. FLAT ROBERTSON- WALKER BACKGROUND 

We first consider the case of a flat RW background, whose metric in conformal coordinates reads | 

ds 2 =a 2 (r/)(-dr/ 2 +dx 2 ) . (23) 

Under the assumption that the graviton quantum vacuum preserves the symmetries of the background metric, we can 
conclude that the tensor M^ vaT has the same symmetries. For RW spacetimes the metric is invariant under spatial 
reflections (due to its homogeneity and isotropy), so that M°'^ k = 0. For the remaining two non- vanishing sets of 
components of the tensor, M * ^ and M % i kl , we use the invariance of the metric under spatial rotations. This implies 
that they can be written in the form 

M ^' =fi(v)9 ij , (24a) 
M ijkl = h( V )(g ik gi l -g u gi k ) , (24b) 

where f\{rf) and /^(f?) are functions of time. Note that the non- vanishing components of the tensor have the same 
form as the components of the Riemann tensor in RW spacetimes, apart from the global factors /i and fi. To 
determine these two functions, we recall that the tensor M^ V(JT is proportional to the two-point function of gravitons, 
which in RW backgrounds can be expressed in terms of the two-point function of a massless minimally coupled 
scalar field ptlBjJ. Therefore the functions fi(ji) and /^(f?) must be proportional to {4> 2 (if)). In order to calculate the 



3 For simplicity we have considered the case of a flat RW background. Our results can be easily generalized to the closed and 
open RW spacetimes. 
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constants of proportionality exactly, we need to go beyond symmetry arguments and we must face the exact evaluation 
of the graviton two-point function G^ vp(J {x, x') of Eq. ( |Tc| ) in a RW gravitational background. The result can be found 
in the literature (see for example Ref. p0[ ). We only need to quote the final result of that reference. In the traceless 
{hP = 0), transverse (V ^/V 1 " = 0), and synchronous gauge (hence h ^ 1 — 0), and assuming the graviton vacuum 
state to be homogeneous and isotropic, the coincidence limit of the two-point function has only spatial non- vanishing 
components, and reads [p0[ 

G l3kl (x,x) = 2^{m l m 3 m k *m l * + m"m 3 *m k m l )\F(x,k)\ 2 . (25) 

k 

The sum over k denotes a sum over a three dimensional set of spatial wave vectors. The complex (spacelike) vector 
77i l (k) is defined as m*(k) = (l/v / 2)[e l 1 (k) + ie|(k)]. The vectors ei(k) and e2(k) are spacelike vectors, such that the 
set {ei(k), e2(k), k} forms an orthonormal basis in the three dimensional spacelike hypersections (here k = k/k and 
k = (k-k) 1 / 2 ). The mode function F is given by F(x,k) = F(r?,x,k) = f(r),k)e ik - x /y/32n 3 kV, where V is a constant 
comoving volume [J. The mode functions F(x, k) and /(k, 77) are a solution to the equations 

□F(x,k) = 0, (26a) 

f( V , k) + —f( v , k) + k 2 f( v , k) = , (26b) 
a 

respectively, i.e., correspond to the dynamical equation of a massless minimally coupled scalar field in a RW back- 
ground. Here the dot denotes derivation with respect to the conformal time variable 77. 

The two-point function, given in Eq. (p5|), can be simplified by making use of the identity e\e 3 l +e l 2 e 3 2 + k~ 2 k l k 3 — g 13 . 
We can then perform the sum over momenta as \F(x, k)| 2 depends only on the modulus of k. The final result reads 

G i3kl (x, x) = 1 Q<fV + \ftf h ~ ffV') (<P(v)) , (27) 

where (4> 2 (ri)) = Ylu. k)| 2 is the coincidence limit of the two-point function of a massless minimally coupled 
scalar field in a RW bac kgrou nd. We c an now insert this formula in Eq. (ph in order to read off the expressions for 



fi{rf) and /^(f?) in Eqs. (24a) and (241). This procedure yields 



f2(v) = ~{4> 2 (v))- (28) 

The quantum correction to the classical EM action due to the coupling with the graviton degrees of freedom can now 
be obtained from Eq. (||) by making use of the above results. We get 

(4m) = / d^v 7 ^ [-4a 2 (r,)/ 1 ( ?7 )^F M + 2f 2 (r,)F ij F i3 ] . (29) 

Note that the one- loop effective action for the electromagnetic field is <Sem + (^em)' wn ich is a divergent quantity 
and has to be suitably renormalized. This is accomplished by the renormalization of (</> 2 (77)), for example, by means 



of adiabatic regularization |22 



Having the EM effective action and assuming that one first solves the pure gravity part in order to get the corrected 
background metric after graviton back reaction (that is, we assume one solves the pure gravitational part and gets 
the new scale factor a(r))), one can get the corrected Maxwell's equations from the variation of the effective action, 
namely, 5S e s/5A^ = 0. We will assume that fi(rj) and ^(f?) have a time variation much slower than that associated 



4 The choice of vacuum corresponds to that used in Ref. |^TJ. This graviton vacuum is homogeneous, isotropic, and the same 
for the two helicity states of the gravitons (+2, —2). 

5 The normalization of the mode function F differs from that of Allen J2FJ] in a k _1 factor. The reason is that we have defined 
the graviton fluctuations via <7 M „ — » <7 M „ + nh^v and our graviton two-point function, (see Eq. (]ic|)), is given in terms of this 
ft M ", whereas Allen has <7 M „ — > g M „ + /i M „, and defines the graviton two-point function in terms of /t M „. 
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with the EM field, so that we can approximate the equations of motion as in Eq. ( fl5| ) . The source-free equations are 
the usual ones 

' V • B = , (30a) 



a{rf) 



a(rj) a(n) a{-q) 

where the dot denotes d/drj, and V is vector notation for d/dx. 
The other two equations read 



1 b + H^) B = * VxE, (30b) 



' e e ff(r?)V • E = , (31a) 



a{rj) 



-L<*ff fa)E + ^e eff („)E = 4t^(»?)V x B . (31b) 
«(»?) «W) «w) 

In the first term on the left hand side of the last equation we have discarded a contribution proportional to the time 
derivative of the effective electric permeability e c s(rf) since, as already discussed in Eq. (|l5|), we are assuming that it 
does not change significantly over the wavelength of the EM field. The effective electric permeability and magnetic 
permitivity tensors in RW backgrounds are proportional to the identity (3 x 3) matrix, namely, e l j = e c g(r))5 l j and 
l^j = fiesWPj, with 

6 cff (r7) = l + 8o 2 (T7)/i(t/) , (32a) 
Atoff (n) = 1 + 8/2(77) . (32b) 

Hence the presence of gravitons introduces a time dependent effective refraction index n e s On) = \J e ff('7)A 1 off (jj) for 
a travelling EM wave, and therefore implies a time dependent speed of light in the medium H 

v cS (v) = 1 - 4[a 2 (r?)/ 1 (r?) + / 2 (r?)] = 1 + ^(<f> 2 ( V )) ■ (33) 

6 

This effective speed of light is the same for all directions of propagation and for all polarisations of the EM radiation 
field, in agreement with the isotropy and homogeneity of RW spacetimes. In Appendix ^ we give an alternative 
derivation of this result based on a direct study of the dispersion relation for light in the graviton modified medium. 

The renormalized two-point function {(j) 2 {rf)) does not have a definite sign (see, for example, Refs. [ p7p^ ]). The 
effective speed of light in the graviton vacuum can be greater or smaller than that in free space, depending on the 
particular form of the scaling parameter a{rf). In any case, the correction is extremely small, typically proportional 
to the ratio of the spacetime scalar curvature and Planck's curvature R/Rp. As we go back in time towards the 
Big Bang singularity, the modulus of the correction to the phase velocity increases. Of course, we cannot trust this 
calculation for such early times since the correction would become too large, and since General Relativity would not 
be valid as an effective low energy/large distance theory in that regime. 

It is worth mentioning that similar results are obtained due to QED vacuum polarisation || . The QED effects are 
generically much larger than the graviton corrections. However, there are situations in which the QED correction 
vanishes, while the graviton correction does not. To show an explicit example, let us consider de Sitter spacetime. 
Virtual e~e + pairs modify Maxwell's equations as follows S 



7aR 



10807TTO 2 



= , (34) 



where m is the mass of the electron and a the fine structure constant. The corrected equations coincide with Maxwell's 
equations up to a trivial normalization and the dispersion relation is the classical one. However, graviton vacuum 
corrections in de Sitter spacetime do affect the propagation of individual light rays. 



This velocity corresponds to the phase velocity v p hase = c/n e g- 
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B. STATIC AND SPHERICALLY SYMMETRIC BACKGROUNDS 



In this section we consider a static and spherically symmetric spacetime described by the metric 



d,s 2 = 



-A{r)dt 2 + B(r)dr 2 + r 2 dQ 2 



(35) 



In order to compute the tensor M^ vaT one needs to calculate the graviton two-point function, evaluated in an 
arbitrary quantum state. It would be a rather formidable task to explicitly calculate such an object. Instead, we will 
use symmetry arguments and assume that the graviton quantum state preserves the symmetries of the background 
metric to derive the basic structure of the tensor M ta " TT . It is shown in Appendix |b] that the tensor M >lvaT can be 
written as 



M ^ aT (r) = f 1 (r)W xv U' TT + h{r)V^V aT + f 3 (r)(X^X aT + Y^Y aT ) + f 4 (r){W^"W (TT + Z» v Z aT ) 



(36) 



where the functions fi(r) depend on the particular choice of vacuum. The antisymmetric order two tensors 
Ui lv ,V' i ' / ,W >tu ,Xi iU ,Y' JlV , and Z» v are defined in Appendix || Just as in the case of RW backgrounds, the only 
non- vanishing components are M 010 - 7 and M yfei , with i,j,k,l spatial indices. In general the structure of the tensor 
M lll " yp in Eq. ([36|) is much more complicated than that for the Riemann tensor corresponding to the metric d35|). 
However, this form for the tensor M pvaT is good enough to carry an analysis parallel to that of Drummond and 
Hathrell §. 

Starting from Eq. ( |l9| ) , which describes the propagation of an EM wave in the presence of gravitons, we show in 
Appendix |b] that it has non-trivial solutions only when the wave momentum satisfies the following determinantal 
condition 



k 2 [(1 + Sf 3 )k 2 - 8(/ 3 + h)k 2 r - 8(/ 2 + h)(k 2 + k 2 )} 

x [(1 + 8/ 3 )(l + 8fr)k 2 - 8(/ 3 + / 4 )(1 + 8A)fc r 2 - 8(/i + / 4 )(1 



0(*? + ^)]=o. 



(37) 



The solution k 2 = corresponds to the usual dispersion relation, in which the light ray follows the null geodesies of the 
background metric. Apart from this (trivial) case, the previous equation admits new dispersion relations, depending 
on the direction of propagation and polarisation of the EM radiation field. 

When the light ray moves radially (kg — k^ — 0) the determinantal condition, (see Eq. (|37j)), has two possible 
solutions 



(1 + 8/ 3 ) (-k 2 + k 2 ,) - 8(/ 3 + h)k 2 = , 
(1 + 8/ 3 )(l + SAX-fc 2 + k 2 ) - 8(/ 3 + J 4 )(l + 8.A)fc 2 = . 



(38a) 
(38b) 



If we assume that (1 + 8/1) ^ the two dispersion relations that follow from the above equations are the same, which 
agrees with the fact that as the gravitational background is spherically symmetric, a radial EM wave should not be 
affected by birefringence. We can conclude that for radial light rays the absolute value of the quantum corrected 
velocity is given by 



h 

When the EM wave moves transversally (k r — kg 

-(l + Sf 3 )k 2 
-(1 + 8/0*1 



1 - 4(/ 4 + h) ■ (39) 

= 0), we get the following two possible solutions (for 1 + 8/3 7^ 0) 

(1 - 8f 2 )k 2 = , (40a) 
(1 - 8.U)k 2 = . (40b) 



As opposed to the previous case, the two dispersion relations that follow from these two equations are different. Light 
rays propagate with different velocities 



1 



4(/ 2 + / 3 ) , 
4(/i + h) , 



(41a) 
(41b) 



depending on their polarisation. Similarly, given the symmetry under the exchange of <f> to 0, azimuthal moving EM 
waves (k r = k<f, = 0) have the same dispersion relations as transverse light rays. Both for traverse and azimuthal 
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propagation we obtain gravitational birefringence due to graviton vacuum fluctuations. Similar results are obtained 
from corrections due to QED vacuum effects. 

It is worth mentioning that the previous results are applicable to any static and spherically symmetric back- 
ground. We will now consider two particular examples, the Schwarzschild background and the Reissner-Nordstrom 
background The Schwarzschild spacetime is described by the metric 

ds 2 = - (l - dt 2 1 dr 2 + r 2 (d6 2 + sin 2 Odcf) . (42) 

In this case, using dimensional analysis, we can write the functions fi(r) (i = 1, . . . , 4) as fi(r) = (Mp / M) 2 Ti(r /2M G) 
with Ti four dimensionless functions and Mp the Planck mass. As we have already mentioned, the exact form of 
these functions is unknown, since it is not possible to calculate explicitly the graviton two-point function for the 
Schwarzschild background. However, near the horizon r « 2MG, we can approximate 

M^ OT (r « 2MG) m (¥f \ [W)U^U aT +r a (l)V> w V aT + T 3 (l)(X^ X aT + Y^Y aT ) 

+ F A (l){W^W aT + Z^Z aT )] . (43) 

Hence, near the horizon the birefringence effects induced by gravitons are of order 0(Mp/M 2 ). 
The Reissner-Nordstrom spacetime describes the metric of a static and charged black hole 

l2 / 2MG g 2 G\ , , / 2MG Q 2 G\~ 1 , 2 2<rwi2 . ,„ lj2x 

ds 2 = - 1 + ^—^ )dt 2 + 1 + 5» dr 2 + r 2 (d6» 2 + sin 2 (9d</> 2 ) . (44) 

\ r 47rr z / \ r 47rr z / 

This spacetime has in general two event horizons r±. The exterior one at r + = MG[1 + (1 — Q 2 /47rAf 2 ) 1 / 2 ] coincides 
with that of the Schwarzschild metric in the limit Q 2 <C M 2 . We can carry out an analysis similar to the previous 
one, and write fi(r) = (Mp/M) 2 Qi(r/2MG,r/QyfG) with Qi four new dimensionless functions. Near the exterior 
event horizon (r = r+) we can approximate 



ir ,r (r«r + )« 



( M P 
V M 



2MG , Q^GJ \2MG' Qy/G 

I (x^x ffT + Y^Y aT ) + g 4 ( -^±-, ^4 



(45) 



Hence, near the outer event horizon the birefringence effects induced by one-loop gravitons are again of order 
0(M 2 /M 2 ). 



V. CONCLUSIONS 



In this paper we have computed the quantum corrections to (classical) Maxwell's equations due to the interaction 
between the EM field and the fluctuations of the spacetime metric. The modified equations look like Maxwell's 
equations in the presence of a linear medium, with electric permitivity and magnetic permeability tensors proportional 
to the coincidence limit of the graviton two-point function. From the modified equations we have found the quantum 
corrections to the dispersion relations. In general, as for the case of linear media in classical electrodynamics, the 
speed of light depends on the direction of propagation and on the state of polarisation. The quantum corrections we 
have computed should be considered along with the graviton back reaction on the background metric (this is crucial, 
for example, for the cancellation of the gauge-breaking dependence; see Refs. |l3|,|l7]] ) . Throughout the paper we have 
assumed that the metric does contain such (back reaction) corrections, and have focused on the interaction of the 
gravitons with the electromagnetic field. 



7 As we already mentioned at the end of Section ||, the background metric should be a solution of the semi-classical Einstein 
equation. This is not the case for the particular examples we will be considering in this Section. We include them only for 
illustrative purposes. 
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We have shown that, when the EM field wavelength A is small compared to the typical scale L c of variation of the 
permitivity and permeability tensors, the corrections to the speed of light are proportional to (Lp / L c ) 2 . We have 
described in detail the quantum corrections in both RW gravitational backgrounds and static spherically symmetric 
spacetimes. 

For RW spacetimes we have computed the quantum corrections by two different methods: the analysis of the 
modified Maxwell's equations in a coordinate basis and the study of the dispersion relations in a local Lorentz frame. 
The corrections we found are similar to those of Ref. Jl?]] , where the analysis was based on the effect of gravitons on 
the spacetime null geodesies. The results agree qualitatively but not quantitatively (as is to be expected), since the 
coupling of gravitons to point particles is in general different to the coupling to massless fields. 

For spherically symmetric spacetimes we have been able to estimate the quantum corrections using symmetry and 
dimensional arguments (see Appendix |b|), avoiding the explicit computation of the graviton two-point function. 

The quantum corrections for small wavelengths are also qualitatively similar to those produced by virtual e _ e + in 
non-trivial backgrounds. However, for some particular cases (as de Sitter background), the QED vacuum does not 
affect the propagation of photons, whereas the graviton vacuum does induce a modification on the propagation of 
light rays. 

The opposite limit, A 3> L c , is more interesting from a phenomenological point of view. Assuming a non-trivial 
and random spacetime structure at scales of order L c , the modified field equations are similar to the ones describing 
the propagation of classical waves in random media. To lowest order, it is possible to describe the wave propagation 
in terms of an "effective medium" . The average corrections to the speed of light are independent of the wavelength 
and proportional to (Lp/L c ) 2 . The "effective medium" is only an averaged description. Stochastic fluctuations are 
expected to occur, for example, in the arrival time of photons coming from point sources |24|,^5]]. To next order, the 
corrections are proportional to L 2 p E /(L C A), where E is the photon energy . In the extreme case Lp ~ L c (which we 
cannot reach within our effective field theory approach) the corrections would be proportional to E/Ep, where Ep 
is the Planck energy. In this regime quantum gravity effects increase with energy, while other medium effects should 
decrease with energy. They can be distinguished by this property and could be relevant in cosmological situations Q . 

In the effective field theory approach to quantum gravity one also expects classical, stochastic fluctuations of the 
spacetime metric. Its dynamics should be described by a "semi-classical Einstein-Langevin equation" ^(|. These 
classical fluctuations will also affect the propagation of photons (see for example |^7]j). 
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APPENDIX A: STRUCTURE OF THE TENSOR M FOR HOMOGENEOUS AND ISOTROPIC 

BACKGROUNDS 

We consider a flat homogeneous and isotropic spacetime described by the metric (conformal coordinates) 

ds 2 = a 2 (? ? )(-d?7 2 + dx 2 + dy 2 + dz 2 ) . (Al) 

In order to compute the tensor M tJ,vcrT one needs to calculate the graviton two-point function, evaluated in an arbitrary 
quantum state. We will assume that the graviton state preserves the symmetries of the background metric. The aim 
of this appendix is to use those symmetries to find the structure of the tensor M^ vrjT given in Eq. (|9|). In order to 
obtain this structure we choose an orthonormal basis of the spacetime under consideration. The orthonormal basis is 
defined by the vector fields 



1 d 

n a(?y) drj 
1 d 
x a(r/) dx 

- 1 9 
Jy a(rj) dy 



(A2a) 
(A2b) 
(A2c) 
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- 1 9 

z a(rj) dz 

We also introduce the following set of antisymmetric tensors 

= e^e" - e u e^ 

w x c r) c a; c r; c a; ) 

u y G T) c y K ri c y > 

TJV" = e^e" - e'V 

^ z c 7; c 2 c »7 c z ) 

x H Z J/ X ' 



(A2d) 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 



Given the fact that the e's form an orthonormal basis, this set of six tensors constitutes a basis for the antisymmetric 
order two tensors. 

It is easy to see that the symmetries of the tensor M^ vaT (antisymmetric in the first two indices, and the last two, 
and symmetric under the exchange of the two pair of indices) imply that the tensor must be a linear combination of 
the following kind (with coefficients that may depend only on r/) 



+ an(vWVr + VrU° T ) + a^U^U^ + a 13 ( V )(U^Vr + V^T) 



y 1 y y > 



a u ( V )(UrV° T + Vy^UD + a 15 (r,)(Urvr + VJTU?) + "16 



y y 

, , .„ , , )+ ai8 (v)(vrvz aT + vrvr) 



axr(v)(vrv° r +vrvr 



(A4) 



Since the metric is homogeneous and isotropic the functions a2, 03, 04, 05, 06, as, ag, aio, an, ai3, ai4, ai5, an, a-is, and 
020 must vanish identically. Here we should stress once again the fact that the quantum state of the gravitons does 
not break these symmetries. The background symmetries also imply that a\ = aj = a\2 = a and ai6 = aig = a2i = (3. 

After these symmetry considerations, we can finally write the general expression of the tensor M A " y<TT for a homo- 
geneous and isotropic background 



(A5) 



This is all the information we can obtain regarding the tensor structure of A4 fJ/l " 7T by making use of the symmetries 
of the gravitational background. The functions a(rj) and (3(r)) will depend on the particular choice of vacuum and are 
not known a priori. 

We now make use of the newly obtained form for M^ var , (see equation (A5)) to solve the equation of motion for 
the EM field in the presence of one-loop quantum fluctuations in flat homogeneous and isotropic backgrounds (see 
Eq. (|l9|)). Following the very same steps as described in Ref. ||, we introduce the six (dependent) functions 



u x fnv^x 
u y = ffivUy 



X = f^V^x > 

'y = fixuVy 1 
'z = ffiuV z ■ 



From Eq. (pi) it follows that 



fliv 



k^(x v 



(A6) 



(A7) 



for some gauge potential a M . Hence, for a given EM wave momentum f^ u has three independent components (one 
amplitude and two polarisations), as we still have the choice of gauge for the EM field. Since is gauge invariant, 
without loss of generality we consider the Coulomb gauge and choose a^ = 0. With this choice the three non- vanishing 
components for the gauge potential are a x , a y , and a 2 , so that 



fr/x 


krj&x 


fxy 


— k x a,y 


kyQ>xi 


fvv 


— k'qCLy 


fyz 


— ky(X Z 


k z dy. 


ff]Z 


= k^a z 


fzx 


= k z a x 


— k x a z , 



(A8) 
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and we can then write 







V x 


— %fyz\ 


Uy 




Vy 


— ^fzx\ 


u z 




v z 





so that in terms of the independent set {u x , u y , u z } the three dependent ones can be written as 

v x = -r~(kyU z k z lly) , 
v y = ~r~i K zU x k x ii z ) , 

V Z = ~j—(k X Uy kyUx) . 

Kyi 



(A9) 

(AlOa) 
(AlOb) 
(AlOc) 



With all these definitions in mind we project equation ( |l9| ) onto the three tensors U x , U y , and U z (that yield u x ,u y , 
and u z , respectively), to obtain the following set of equations 



= k 2 u x + 8a[u x (-k 2 + kl) + k x (k y u y + k z u z )] + 8/3k v (k y v z - k z v y ) , 
= k 2 u y + 8a[u y (-k 2 + k 2 y ) + k y (k x u x + k z u z )] + 8(3k v (k z v x - k x v z ) , 
= k 2 u z + 8a[u z (—k 2 + k z ) + k z (k x ii x + k y u y )] + 8(3k v (k x v y — k y v x ) . 



(Alia) 
(Allb) 
(Allc) 



We point out that the components of the vector k correspond to the orthonormal basis given in Eqs. (A2), so that 
k = k n e v + k x e x + k y e y + k z e z . If we write the components (v x , v y ,v z ) as functions of (u x , u y ,u z ), we have 



= k 2 u x + 8a[u x (-k 2 + kl) + k x {k y u y + k z u z )] + 8/3[k y (k x u y - k y u x ) - k z (k z u x - k x u z )) , 

= k 2 Uy + 80t[U y (-k 2 + ky) + ky{k X U X + k Z U Z )] + 8/3[k Z (kyU Z - k Z Uy) - k X (k X Uy ~ kyU X )\ , 

= k 2 u z + 8a[u z (-k 2 + k 2 z ) + k z (k x u x + k y u y )] + 8j3[k x {k z u x - k x u z ) - k y (k y u z - k z u y )] . 



(A12a) 
(A12b) 
(A12c) 



Imposing the condition that the determinant of this set of equations vanishes, so that we do not obtain the trivial 
solution (u x = Uy = u z = 0), we obtain the following determinantal restriction 

fc 2 (l + 8a) [(1 + 8a)k 2 - 8(a + (3)(k 2 x + k 2 y + k 2 z )] 2 = . (A13) 

Let us assume that the EM radiation is characterized by a three-dimensional momentum k, such that k 2 + k 2 + k 2 = 
k • k. The non-trivial dispersion relation becomes then 



(1 + 8a) (-^ + k ■ k) - 8(a + /3)k ■ k = . 
The EM radiation waves will propagate with the following dispersion relation 



s = 1 



k k 



(a + /3), 



or equivalently 



1 - 4(a + 0) 



(A14) 



(A15) 



(A16) 



We can now comp are with the previous formulation of this problem in terms of the two-point function of the gravitons 
(see section IV A). We already know from symmetry considerations that the tensor M^ vaT can be written as 



M^ aT = a(r])(U^U^ T + UfU" + U^U° T ) + (3(r 1 ){V^V x aT + V* v V? r + V z ^V z aT ) . 



(A17) 



Let us now calculate the only non- vanishing components of this tensor (M ??l ' u and M l ^ mn ) in the coordinate basis 
defined by the vector fields d v ,d x , d y and d z . We get 
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M ijmn = p(ri){V* j V™ n + V^V™ n + V l JV™ n ) = (i{vW m g ]n ~ g m 9 3m ) , 



(A18) 
(A19) 



which means in particular that the functions a(rj) and P(rj) introduced in this appendix are related to the functions 
defined in Eqs. (24a) and (24b) as fi(r]) — a(r}) / a 2 (rf) and f 2 (r/) = Hence, we conclude that 



= 1 - 4(a + 0) = 1 - i[a 2 ( V )fi(v) + f 2 (v)} = 1 + y (<?(j,)) > 



(A20) 



which agrees with our previous result obtained from Maxwell's equations and the two-point function of the gravitons 
(see equation (R3I) ) . 



APPENDIX B: STRUCTURE OF THE TENSOR M FOR STATIC AND SPHERICALLY SYMMETRIC 

BACKGROUNDS 



We assume a static and spherically symmetric spacetime described by the metric 



ds 2 = ~A{r)dt 2 + B(r)dr 2 



(Bl) 



In order to compute the general form of the tensor M Miy<TT we will follow the same steps as in Appendix |A|. We choose 
an orthonormal basis of the spacetime given by the vectors e t , e r , eg, and e^. This orthonormal basis is defined by 



e* = 



d 



Or 



ld_ 

rd9 ' 
1 d 



r sin 6 dcj) 

We also introduce the following basis of antisymmetric order two tensors 



(B2a) 
(B2b) 
(B2c) 
(B2d) 





— *=* e r 


_ e *V< 






v u 




— e i e 6» 


- e v e M 






v u 




— p*V^ 

— e r e e 






- e^e 1 ! 


- P y P M 



(B3a) 
(B3b) 
(B3c) 
(B3d) 
(B3e) 
(B3f) 

The tensor M^ vaT must be a linear combination of the following kind (with coefficients that may depend only on r) 

M^ aT = ai {r)U^U aT + a 2 {r){U^V aT + V^U aT ) + a 3 {r){U^ X aT + X^U aT ) 

+ a 4 (r)(U^Y aT + Y^U ar ) + a 5 (r){U^W' TT + W^U^) + a 6 (r)(U^ Z aT + Z^U aT ) 

+ a 7 (r)V^V lTT + a 8 {r){V^X aT + X» V V T ) + a^r^V^Y^ + Y^V aT ) 

+ a X0 {r){V tlv W' TT + W^V aT ) + ail {r){V^ Z aT + Z^V aT ) + a l2 {r)X» v X aT 

+ a^r^X^Y^ + Y^X aT ) + ai 4 (r)(X^W ffT + X aT ) + a 15 (r){X^ Z aT + Z^X aT ) 

+ a w (r)Y^Y aT + a l7 {r){Y^W aT + W^Y aT ) + a ls (r)(Y^Z aT + Z^Y aT ) 

+ a w (r)W^W aT + a 2a {r){W^Z aT + Z^W aT ) + a 21 (r)Z^ Z aT . (B4) 



Since the metric is static, we can make use of the time inversion invariance to show that the terms in Eq. (B4) 
involving the functions a 2 , a§, ag, as, ag, an, ai5, an, and ai$ must vanish identically Here we should stress once 
again that the quantum state of the gravitons does not break time inversion invariance. (This is not true in general. 
For example, the Unruh vacuum state in Schwarzschild spacetime does break the time inversion symmetry.) 



13 



Because of the spherical symmetry of the metric (spatial inversion as well) the following coefficients have to vanish: 
a 3 , 0,4, dio, on, ^13, and a 2 Q. Furthermore, the coefficients an and die must be equal as they have to be invariant 
under spatial rotations. The same is true for the pair a±g and 021- 

After these symmetry considerations, we can finally write the general expression for the tensor M'"" TT in a static, 
spherically symmetric background 



M"" ffI = /i(r)C/^C/ CTT + f 2 {r)V tlv V' JT + f 3 (r){X" p ' X aT + Y" u Y aT ) + f A (r){W' i "W' Tr + Z^Z aT ) 



(B5) 



The functions fi(r) (with i — 1, ... ,4) will depend on the particular choice of graviton state and are not known a 
priori. 

As in Apendix A we introduce the six functions 



x = f^X^ 



(B6) 



In the Coulomb gauge, the three non-vanishing components for the gauge potential are a r , ag, and a^, so that 



and we can then write 



ftr — fc(O r 

ft<p = kta<f, 



U = 2 f tr 

y = Vt<t> 

z — 2f r( p 



fte = k t ae; 
f r g = k r a e - k e a 
TBih = kgad, - 



x = 2 f te ; 
w = 2f r g; 
v -- 



(B7) 



so that in terms of the independent set {it, x, y}, the three dependent ones can be written as 

w = —{k r x — kgu) , 

h 

z = -r(k r y - k^u) , 
v = -}-(kgy - k^x) . 



(B8) 

(B9a) 
(B9b) 
(B9c) 



We now project equation ( ^9[ ) onto the three tensors U, X, and Y (that yield u, x, and y, respectively), to obtain the 
following set of equations 



= k 2 u - 8fiu(k 2 - k 2 ) + 8f 3 k r (xkg + ykj,) + 8fik t {wkg + zk^) , 
= k 2 x + 8fiuk r kg + 8f 2 vk t k (f> - 8f 3 x(k 2 - kg) + 8f 3 ykgk (t> - 8f4wk t k r , 
= k 2 y + 8fiuk r k$ - 8f 2 vk t kg + 8f 3 xkgk tj) - 8f 3 y(k 2 - kfy - 8f±zk t k r . 

If we write the components (v, w, z) as functions of (it, x, y), we have 

= [k 2 (l + 8/1) - 8(/i + f 4 )(k 2 + k 2 )} u + 8(/ 3 + U)k r k e x + 8(f 3 + U)k r k^y , 
- 8(/i + h)k r k e u + [k 2 (l + 8/3) - 8(/ a + f 3 )kl - 8(f 3 + U)k 2 ] x + 8(f 3 + U)kgk^y , 
- 8(/i + h)k r k^u + 8(/ a + f 3 )k e k^x + [k 2 (l + 8/3) - 8(/ 2 + f 3 )k 2 - 8(/ 3 + U)k 2 } y. 

In this case the components of the vector k in the orthonormal basis are given by k = fc'ej + k r e r + k e eg + k^e^. 
Imposing that the determinant of the set of three equations vanishes, we obtain the following detcrminantal condition 

k 2 [(1 + 8/ 3 )fc 2 - 8(/ 3 + h)k 2 - 8(f 2 + f 3 ){k 2 + k 2 )] 

x [(1 + 8/ 3 )(l + 8h)k 2 - 8(/ 3 + / 4 )(1 + 8h)k 2 - 8(/i + U)(l + 8f 3 )(k 2 + k 2 )} = 0. (B12) 



(BlOa) 
(BlOb) 
(BlOc) 

(Blla) 
(Bllb) 
(Bile) 
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